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We introduce a parameterized high-density equation of state (EOS) in order to systematize the 
study of constraints placed by astrophysical observations on the nature of neutron-star matter. To 
obtain useful constraints, the number of parameters should be smaller than the number of neutron- 
star properties that have been measured or will have been measured in the next several years. 
And the set must be large enough to accurately approximate the large set of candidate EOSs. 
We find that a parameterized EOS based on piecewise polytropes with 3 free parameters matches 
to about 4% rms error an extensive set of candidate EOSs at densities below the central density 
of 1.4 Mq stars. Adding observations of more massive stars constrains the higher density part of 
the EOS and requires an additional parameter. We obtain constraints on the allowed parameter 
space set by causality and by present and near-future astronomical observations. In particular, 
we emphasize potentially stringent constraints on the EOS parameter space associated with two 
measured properties of a single star; and we find that a measurement of the moment of inertia of 
PSR J0737-3039A can strongly constrain the maximum neutron-star mass. We also present in an 
appendix a more efficient algorithm than has previously been used for finding points of marginal 
stability and the maximum angular velocity of stable stars. 
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I. INTRODUCTION 

Because the temperature of neutron stars is far below 
the Fermi energy of their constituent particles, neutron- 
star matter is accurately described by the one-parameter 
equation of state (EOS) that governs cold matter above 
nuclear density. The uncertainty in that EOS, however, 
is notoriously large, with the pressure p as a function 
of baryon mass density p uncertain by as much as an 
order of magnitude above nuclear density. The phase 
of the matter in the core of a neutron star is similarly 
uncertain: Current candidates for the EOS include non- 
relativistic and relativistic mean-field models; models for 
which neutron- star cores are dominated by nucleons, by 
hyperons, by pion or kaon condensates, and by strange 
quark matter (free up, down, and strange quarks); and 
one cannot yet rule out the possibility that the ground 
state of cold matter at zero pressure might be strange 
quark matter and that the term "neutron star" is a mis- 
nomer for strange quark stars. 

The correspondingly large number of fundamental 
parameters needed to accommodate the models' La- 
grangians has meant that studies of astrophysical con- 
straints (see, for example, [TJ [2 El HI E] and references 
therein) present constraints by dividing the EOS can- 
didates into an allowed and a ruled-out list. A more 
systematic study, in which astrophysical constraints are 
described as constraints on the parameter space of a pa- 
rameterized EOS, can be effective only if the number of 
parameters is smaller than the number of neutron-star 
properties that have been measured or will have been 
measured in the next several years. At the same time, 
the number of parameters must be large enough to accu- 
rately approximate the EOS candidates. 



A principal aim of this paper is to show that, if one uses 
phenomenological rather than fundamental parameters, 
one can obtain a parameterized EOS that meets these 
conditions. We exhibit a parameterized EOS, based on 
specifying the stiffness of the star in three density inter- 
vals, measured by the adiabatic index F — d\og P/dlog p. 
A fourth parameter translates the p(p) curve up or down, 
adding a constant pressure - equivalently fixing the pres- 
sure at the endpoint of the first density interval. Finally, 
the EOS is matched below nuclear density to the (pre- 
sumed known) EOS. An EOS for which F is constant is 
a polytrope, and the parameterized EOS is then piece- 
wise polytropic. A similar piecewise-polytropic EOS was 
previously considered by Vuille and Ipser [6^; and, with 
different motivation, several other authors [3 HI [H ITU] 
have used piecewise polytropes to approximate neutron- 
star EOS candidates. In contrast to this previous work, 
we use a small number of parameters chosen to fit a wide 
variety of fundamental EOSs, and we systematically ex- 
plore a variety of astrophysical constraints. Like most of 
the previous work, we aim to model equations of state 
containing nuclear matter (possibly with various phase 
transitions) rather than pure quark stars, whose EOS is 
predicted to be substantially different. 

As we have noted, enough uncertainty remains in the 
pressure at nuclear density, that one cannot simply match 
to a fiducial pressure at Pnuc- Instead of taking as one 
parameter the pressure at a fiducial density, however, 
one could match to the pressure of the known subnuclear 
EOS at, say, 0.1 p^uc and then use as one parameter a 
value of Fq for the interval between O.lpnuc and pnuc- 
Neutron-star observables are insensitive to the EOS be- 
low pnuc, because the fraction of mass at low density is 
small. But the new parameter Fq would indirectly affect 
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observables by changing the value of the pressure at and 
above nuclear density, for fixed values of the remaining 
r^. By choosing instead the pressure at a fixed density 
Pi > Pnuc, we obtain a parameter more directly con- 
nected to physical observables. In particular, as Lattimer 
and Prakash j5] have pointed out, neutron-star radii are 
closely tied to the pressure somewhat above nuclear den- 
sity, and the choice pi = p{pi) is recommended by that 
relation. 

In general, to specify a piecewise polytropic EOS with 
three density intervals above nuclear density, one needs 
six parameters: two dividing densities, three adiabatic 
indices Ti , and a value of the pressure at an endpoint of 
one of the intervals. Remarkably, however, we find (in 
Sec. [v]) that the error in fitting the collection of EOS 
candidates has a clear minimum for a particular choice 
of dividing densities. With that choice, the parametrized 
EOS has three free parameters, ri,r2 and pi, for densi- 
ties below 10^^ g/ar? (the density range most relevant 
for masses IAMq), and four free parameters (an ad- 
ditional Fs) for densities between 10^^ g/av? and the 
central density of the maximum mass star for each EOS. 

With the parameterization in hand, we examine in 
Sect . |Vl| astrophysical constraints on the parameter space 
beyond the radius-pi relation found by Lattimer and 
Prakash [5] . Our emphasis in this first study is on present 
and very near-future constraints: those associated with 
the largest observed neutron-star mass and spin, with a 
possible observation (as yet unrepeated) of neutron-star 
redshift, with a possible simultaneous measurement of 
mass and radius, and with the expected future measure- 
ment of the moment of inertia of a neutron star with 
known mass. A companion paper will investigate 
constraints obtainable with gravitational-wave observa- 
tions in a few years. The constraints associated with the 
largest observed mass, spin, and redshift have a similar 
form, each restricting the parameter space to one side of 
a surface: For example, if we take the largest observed 
mass (at a 90% confidence level) to be 1.7 Mq, then 
the allowed parameters correspond to EOSs whose max- 
imum mass is at least 1.7 Mq. We can regard M^^^ as 
a function on the 4-dimensional EOS parameter space. 
The subspace of EOSs for which A^fmax = 1.7M0 is then 
described by a 3-dimensional surface, and constraint is 
a restriction to the high-mass side of the surface. Simi- 
larly, the observation of a 716 Hz pulsar restricts the EOS 
parameter space to one side of a surface that describes 
EOSs for which the maximum spin is 716 Hz. Thus we 
can produce model-independent extended versions of the 
multidimensional constraints seen in [12] . 

The potential simultaneous observation of two prop- 
erties of a single neutron star (for example, moment of 
inertia and mass) would yield a significantly stronger con- 
straint: It would restrict the parameter space not to 
one side of a surface but to the surface itself. And a 
subsequent observation of two different parameters for 
a different neutron star would then restrict one to the 
intersection of two surfaces. We exhibit the result of 



simultaneous observations of mass and moment of iner- 
tia (expected within the next decade for one member of 
the binary pulsar J0737-3039 [131 E]) and of mass and 
radius. Gravitational-wave observations of binary inspi- 
ral can again measure two properties of a single star: 
both mass and an observable roughly described as the 
final frequency before plunge (the departure of the wave- 
form from a point-particle inspiral); and related work in 
progress examines the accuracy with which one can ex- 
tract EOS parameters from interferometric observations 
of gravitational waves in the inspiral and merger of a 
binary neutron star system |llj . 

Conventions: We use cgs units, denoting rest-mass 
density by p, and (energy density)/c^ by e. We define 
rest-mass density as p = m-QU where — 1.66 x 10~^^ g 
and n is the baryon number density. In Sec. |III| however, 
we set c = 1 to simplify the equations and add a footnote 
on restoring c. 

II. CANDIDATES 

A test of how well a parametrized EOS can approx- 
imate the true EOS of cold matter at high density is 
how well it approximates candidate EOSs. We consider 
a wide array of candidate EOSs, covering many different 
generation methods and potential species. Because the 
parametrized EOS is intended to distinguish the parts of 
parameter space allowed and ruled out by present and fu- 
ture observations, the collection includes some EOSs that 
no longer satisfy known observational constraints. Many 
of the candidate EOSs were considered in Refs. [5l [T2l[T4] : 
and we call them by the names used in those papers. 

For plain npep nuclear matter, we include: 

• Two potential-method EOSs (PAL6 [TF and 
SLy il6j), 

• eight variational-method EOSs (APl-4 [17], 
EPS [18;, and WFFl-3 [19 ), 

• one nonrelativistic (BBB2 [5U]) and three rela- 
tivistic (BPAL12 EI], ENG [5g and MPAl [23]) 
Brueckner-Hartree-Fock EOSs, and 

• three relativistic mean field theory EOSs (MS 1-2 
and one we call MSlb, which is identical to MSI 
except with a low symmetry energy of 25 MeV j24j ) . 

We also consider models with hyperons, pion and kaon 
condensates, and quarks, and will collectively refer to 
these EOSs as K/Tr/H/q models. 

• one neutron-only EOS with pion condensates 
(PS [25]), 

• two relativistic mean field theory EOSs with kaons 
(GSl-2 [21]), 

• one effective potential EOS including hyperons 
(BGNIHI [27]), 
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• eight relativistic mean field theory EOSs with hy- 
perons (GNH3 [28 and seven variants Hl-7 [I2j . 

• one relativistic mean field theory EOS with hyper- 
ons and quarks (PCL2 [25]), and 

• four hybrid EOSs with mixed APR nuclear mat- 
ter and colour-flavor- locked quark matter (ALFl-4 
with transition density pc and interaction parame- 
ter c given by pc — 2no, c = 0; pc = 3no, c = 0.3; 
Pc = 3no, c = 0.3; and = 4.5no, c = 0.3 respec- 
tively 

The tables are plotted in Fig. [ijto give an idea of the 
range of EOSs considered for this parameterization. 




log( p in g/cnT*) 

FIG. 1: Pressure versus rest mass density for the set of can- 
didate EOS tables considered in the parameterization. 



III. PIECEWISE POLYTROPE 



A polytropic EOS has the form, 
P{P) = Kp^, 



(1) 



with p the rest-mass density and F the adiabatic index, 
and with energy density e fixed by the first law of ther- 
modynamics,^ 



-p d- 



,1 



(2) 



^ In this section, for simplicity of notation, c = 1. To rewrite the 
equations in cgs units, replace p and K in each occurrence by 
p/c^ and Kji? . Both e and p have units g/cm^. 



For p of the form ([IJ , Eq. ([2]) has the immediate integral 



r-i 



(3) 



where a is a constant; and the requirement lim - = 1 

p^a p 

implies a = and the standard relation 



1 



e = p 



F- r 



(4) 



The parameterized EOSs we consider are piecewise 
polytropes above a density po, satisfying Eqs. ([T|) and 
([3]) on a sequence of density intervals, each with its own 
Ki and F^: An EOS is piecewise polytropic for p > po 
if, for a set of dividing densities po < Pi < P2 < • • • , the 
pressure and energy density are everywhere continuous 
and satisfy 



p{p) ^ K,p^\ d-^-pd-, 
P P 

Then, for F =^ 1, 

e(p) = (1-f a,)p + 



Pi-i < P< P; 



F,:-l 



with 



CLi = 1 

Pi-i 



K, 



F,:-l 



P^^l ■ 



(5) 



(6) 



(7) 



The specific enthalpy^ h is defined as {e + p)/p and is 
given in each density interval by 



h{p) = 1 + ai + 
The internal energy e = | — 

e(p) = flj + 
and the sound velocity Vg is 

Vs{p) = 



K,P 



1 is then 

K.. 



F,;-l' 




(8) 



(9) 



(10) 



^ A note on terminology: When the entropy vanishes, the specific 
enthalpy, h = {e+p)/p, and Gibbs free energy, g = (e+p) / p—Ts, 
coincide. For nonzero entropy, it is the term gdMo or, equiva- 
lently, /idN that appears in the first law of thermodynamics, 
where fi = g/m^ is the chemical potential. Because h = {^+p)/p 
always has the meaning of enthalpy, and because for isentropic 
stars and isentropic flows the constancy of injection energy and 
Bernoulli's law, respectively, are commonly stated in terms of 
enthalpy (see, for example 1311 ). we refer to {t + p)/p as the spe- 
cific enthalpy, rather than the specific Gibbs free energy or the 
chemical potential. 
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Each piece of a piece wise poly tropic EOS is specified by 
three parameters: the initial density, the coefficient Ki, 
and the adiabatic index r^. However, when the EOS at 
lower density has already been specified up to the chosen 
Pi, continuity of pressure restricts Ki+i to the value 



p{Pi) 



(11) 



Thus each additional region requires only two additional 
parameters, pi and Ti+i. Furthermore, if the initial den- 
sity of an interval is chosen to be a fixed value for the 
parameterization, specifying the EOS on the density in- 
terval requires only a single additional parameter. 



IV. FITTING METHODS 

As noted above, in choosing a parameterization for the 
EOS space, our goal was to maintain high enough resolu- 
tion, with a small number of parameters, that a point in 
parameter space can accurately characterize the EOS of 
neutron-star matter. A measure of this resolution is the 
precision with which a point in parameter space can fit 
the available candidate EOSs. We describe in this section 
how that measure is defined and computed for a choice 
of parameter space - that is, for a choice of the set of free 
parameters used to specify piecewise polytropes. 

There is general agreement on the low-density EOS 
for cold matter, and we adopt the version (SLy) given 
by Douchin and Haensel (TB]. Substituting an alter- 
native low-density EOS from, for example, Negele and 
Vautherin [32], alters by only a few percent the observ- 
ables we consider in examining astrophysical constraints, 
both because of the rough agreement among the candi- 
date EOSs and because the low density crust contributes 
little to the mass, moment of inertia, or radius of the 
star. 

Each choice of a piecewise polytropic EOS above nu- 
clear density is matched to this low-density EOS. The 
way in which the match is done is arbitrary, and, again, 
the small contribution of the low-density crust to stellar 
observables means that the choice of match does not sig- 
nificantly alter the relation between astrophysical observ- 
ables and the EOS above nuclear density. In our choice 
of matching method, the first (lowest-density) piece of 
the piecewise polytropic p{p) curve is extended to lower 
densities until it intersects the low-density EOS, and the 
low-density EOS is used at densities below the intersec- 
tion point. This matching method has the virtue of pro- 
viding monotonically increasing EOSs p = p{p) without 
introducing additional parameters. The method accom- 
modates a region of parameter space larger than that 
spanned by the collection of candidate EOSs. It does, 
however, omit EOSs with values of pi and Fi that are 
incompatible, for which the slope of the \ogp vs logp 
curve is too shallow to reach the pressure pi from the 
low-density part of the EOS. 



The accuracy with which a piecewise polytrope 
{pi,Ki,Ti}, approximates a candidate EOS is measured 
by the root mean square or rms residual of the fit to m 
tabulated points Pj,Pj'- 



\ 



( 



\ 



E E (logpj - log + F, logpjf 

\ i 3 

\ Pi<Pi<Pi+l 



(12) 



In each case, we compute the residual only up to the max- 
imum density Pmax that can occur in a stable star - the 
central density of the maximum mass nonrotating model 
based on the candidate EOS. Because astrophysical ob- 
servations can, in principle, depend on the high-density 
EOS only up to the value of Pmax for that EOS, only the 
accuracy of the fit below pmax for each candidate EOS is 
relevant to our choice of parameter space. 

For a given parameterization, we find for each can- 
didate EOS the smallest value of the residual over the 
corresponding parameter space and the parameter val- 
ues for which it is a minimum. In particular, the MIN- 
PACK nonlinear least-squares routine LMDIF, based on 
the Fevenberg-Marquardt algorithm, is used to minimize 
the sum of squares of the difference between the log- 
arithm of the pressure-density points in the specified 
density range and the logarithm of the piecewise poly- 
trope formula, which is a linear fit in each region to the 
log(p)[log/9] curve of the candidate EOS. The nonlinear 
routine allows the dividing points between regions to be 
varied. 

Even with a robust algorithm, the nonlinear fitting 
with varying dividing densities is sensitive to initial con- 
ditions. Multiple initial parameters for free fits are con- 
structed using fixed- region fits of several possible dividing 
densities, and the global minimum of the resulting resid- 
uals is taken to indicate the best fit for the candidate 
EOS. 

We begin with a single polytropic region in the core, 
specified by two parameters: the index Fi and a pressure 
Pi at some fixed density. Here, with a single polytrope, 
the choice of that density is arbitrary; for more than 
one polytropic piece, we will for convenience take that 
density to be the dividing density pi between the first 
two polytropic regions. Changing the value of pi moves 
the polytropic p{p) curve up or down, keeping the loga- 
rithmic slope Fi = dlogp/ dlog p fixed. The low-density 
SFy EOS is fixed, and the density po where the poly- 
tropic EOS intersects SFy changes as pi changes. The 
polytropic index Ki is determined by Eq. (11). This is 
referred to as a one free piece fit. 

We then consider two polytropic regions within the 
core, specified by the four parameters {pi, Fi, pi, F2}, 
as well as three polytropic regions specified by the six 
parameters {pi, Fi, pi, F2, ps, F3}, where, in each case, 
Pi = p{pi). Again changing pi translates the piecewise- 
polytropic EOS of the core up or down, keeping its shape 
fixed. While some EOSs are well approximated by a sin- 
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gle high-density poly trope, others require three pieces to 
capture the behaviour of phase transitions at high den- 
sity. 

The six parameters required to specify three free poly- 
tropic pieces seems to push the bounds of what may 
be reasonably constrained by a small set of astrophys- 
ical measurements. We find, however, that the collection 
of candidate EOSs has a choice of dividing densities for 
which the residuals of the fit exhibit a clear minimum. 
This fact allows us to reduce the number of parameters 
by fixing the densities that delimit the polytropic regions 
of the piecewise polytrope. A three fixed piece fit, using 
three polytropic regions but with fixed pi — 10^^ '' g/cm"^ 
and p2 = lO"'^^ ^ g/cm'^, is specified by the four param- 
eters pi, Fi, and T3. The choice of pi and p2 is 



discussed in Section VB Note that the density of de- 
parture from the fixed low-density EOS is still a fitted 
parameter for this scheme. 



V. BEST FITS TO CANDIDATE EOS 

A. Accuracy of alternative parameterizations 

The accuracy of each of the alternative choices of pa- 
rameters discussed in the last Section, measured by the 
rms residual of Eq. (12), is portrayed in Table |T] The Ta- 
ble lists the average and maximum rms residuals over the 
set of 34 candidate EOSs. 



polytropic pieces models the behaviour of candidates 
well; but this kind of four-parameter EOS does not accu- 
rately fit EOSs with hyperons, kaon or pion condensates, 
and/or quark matter. Many require three polytropic 
pieces to capture the stiffening around nuclear density, 
a subsequent softer phase transition, and then final stiff- 
ening. On the other hand, the six parameters required 
to specify three free polytropic pieces exceeds the bounds 
of what may be reasonably constrained by the small set 
of model-independent astrophysical measurements. As 
mentioned in the introduction, however, an alternative 
four parameter fit can be made to all EOSs if the transi- 
tion densities are held fixed for all candidate EOSs. The 
choice of fixed transition densities, and the advantages of 
such a parameterization over the two free piece fit, are 
discussed in the next subsection. 

The hybrid quark EOS ALF3, which incorporates a 
QCD correction parameter for quark interactions, ex- 
hibits the worst-fit to a one-piece polytropic EOS with 
residual 0.111, to the three-piece fixed region EOS with 
residual 0.042, and to the three-piece varying region 
EOS with residual 0.042. It has a residual from the 
two-piece fit of 0.044, somewhat less than the worst fit 
EOS, BGNIHI, an effective-potential EOS that includes 
all possible hyperons and has a two-piece fit residual of 
0.056. 



B. Fixed region fit 



TABLE I: Average residuals resulting from fitting the set of 
candidate EOSs with various types of piecewise polytropes. 
Free fits allow dividing densities between pieces to vary. The 
fixed three piece fit uses lO^'*'^ g/cm"^ or roughly 1.85pnuc and 
10^^ ° g/cm^ or 3.70p„uc for all EOSs. Tabled are the RMS 
residuals of the best fits averaged over the set of candidates. 
The set of 34 candidates includes 17 candidates containing 
only npep matter and 17 candidates with hyperons, pion or 
kaon condensates, and/or quark matter. Fits are made to 
tabled points in the high density region between 10^* '^ g/cm^ 
or 0.74pnuc and the central density of a maximum mass TOV 
star calculated using that table. 



Type of fit 



All npe^j, K/n/h/q 



Mean RMS residual 



One free piece 0.0386 0, 

Two free pieces 0.0147 0, 

Three fixed pieces 0.0127 0, 

Three free pieces 0.0071 0, 



0285 0.0494 

0086 0.0210 

0098 0.0157 

0056 0.0086 



Standard deviation of RMS residual 



One free piece 0.0213 0, 
Two free pieces 0.0150 0, 
Three fixed pieces 0.0106 0, 
Three free pieces 0.0081 0, 



0161 0.0209 

0060 0.0188 

0063 0.0130 

0039 0.0107 



For nucleon EOSs, the four-parameter fit of two free 



A good fit with a minimal number of parameters is 
found for three regions with a division between the first 
and second pieces fixed at pi = 10^* '' g/cm"^ = l.SSpnuc 
and a division between the second and third pieces fixed 
at p2 = 10^^ ° g/cm^. The EOS is specified by choosing 
the adiabatic indices {Fi,F2,F3} in each region, and the 
pressure pi at the first dividing density, pi = pi{pi). 
A diagram of this parameterization is shown in Fig. [2] 
For this 4-parameter EOS, best fit parameters for each 
candidate EOS give a residual of 0.043 or better, with 
the average residual over 34 candidate EOSs of 0.013. 

The dividing densities for our parameterized EOS were 
chosen by minimizing the rms residuals over the set of 
34 candidate EOSs. For two dividing densities, this is a 
two-dimensional minimization problem, which was solved 
by alternating between minimizing average rms residual 
for upper or lower density while holding the other den- 
sity fixed. The location of the best dividing points is 
fairly robust over the subclasses of EOSs, as illustrated 
in Fig. d 

With the dividing points fixed, taking the pressure pi 
to be the pressure at pi = 1.85pnuc, is indicated by empir- 
ical work of Lattimer and Prakash [5] that finds a strong 
correlation between pressure at fixed density (near this 
value) and the radius of 1.4M0 neutron stars. This choice 
of parameter allows us to examine (in Sec. VIE) the re- 
lation between pi and the radius; and we expect a sim- 
ilar correlation between pi and the frequency at which 
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FIG. 2: The fixed-region fit is parameterized by adiabatic 
indices {Fi , F2 , F3} and by the pressure pi at the first dividing 
density. 
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FIG. 3: Subsets of EOSs with and without kaons, hyperons, 
meson condensates, or quarks, show a fairly robust choice 
of dividing densities whose fit to the candidate EOSs mini- 
mizes residual error. The mean plus one standard deviation of 
residuals for each subset of candidate EOSs is plotted against 
the choice of lower and upper dividing densities pi and p2. 
The left curves show mean residual versus pi with p2 fixed 
at lO^^'" g/cm^. The right three curves show mean residual 
versus p2, with pi fixed at 10^*'^ g/cm^. 



neutron-star inspiral dramatically departs from point- 
particle inspiral for neutron stars near this mass. 

The following considerations dictate our choice of the 



four-parameter space associated with three polytropic 
pieces with two fixed dividing densities. First, as men- 
tioned above, we regard the additional two parameters 
needed for three free pieces as too great a price to pay 
for the moderate increase in accuracy. The comparison, 
then, is between two four-parameter spaces: polytropes 
with two free pieces and polytropes with three pieces and 
fixed dividing densities. 

Here there are two key differences. Observations of 
pulsars that are not accreting indicate masses below 1.45 
Mq (see Sec. VI I, and the central density of these stars 
is below p2 for almost all EOSs. Then only the three 
parameters {pi,ri,r2} of the fixed piece parameteriza- 
tion are required to specify the EOS for moderate mass 
neutron stars. This class of observations can then be 
treated as a set of constraints on a 3-dimensional param- 
eter space. Similarly, because maximum-mass neutron 
stars ordinarily have most matter in regions with densi- 
ties greater than the first dividing density, their structure 
is insensitive to the first adiabatic index. The three piece 
parameterization does a significantly better job above p2 
because phase transitions above that density require a 
third polytropic index T3. If the remaining three param- 
eters can be determined by pulsar observations, then ob- 
servations of more massive, accreting stars can constrain 

The best fit parameter values are shown in Fig. |4] and 
listed in Table |TTT] of Appendix |C] The worst fits of the 
fixed region fit are the hybrid quark EOSs ALFl and 
ALF2, and the hyperon-incorporating EOS BGNIHI. 
For BGNlHl, the relatively large residual is due to the 
fact that the best fit dividing densities of BGNIHI differ 
strongly from the average best dividing densities. Al- 
though BGNIHI is well fit by three pieces with float- 
ing densities, the reduction to a four-parameter fit limits 
the resolution of EOSs with such structure. The hybrid 
quark EOSs, however, have more complex structure that 
is difficult to resolve accurately with a small number of 
polytropic pieces. Still, the best-fit polytrope EOS is able 
to reproduce the neutron star properties predicted by the 
hybrid quark EOS. 

In Appendix [C] Table |III| also shows the various neu- 
tron star structure characteristics for each EOS compared 
to the values of the best-fit piecewise polytrope param- 
eterization for the core. The mean error and standard 
deviation for each characteristic is also listed. 



VI. ASTROPHYSICAL CONSTRAINTS ON 
THE PARAMETER SPACE 



Adopting a parameterized EOS allows one to phrase 
each observational constraint as a restriction to a sub- 
set of the parameter space. In subsections A-D we find 
the constraints imposed by causality, by the maximum 
observed neutron-star mass and the maximum observed 
neutron-star spin, and by a possible observation of gravi- 
tational redshift. We then examine, in subsection E, con- 
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FIG. 4; Parameterized EOS fits to tiie set of 34 candidate 
EOS tables. Tfiere are 17 EOSs witli only ordinary nuclear 
matter (n,p,e,/i); 9 have only hyperons in addition to ordinary 
matter; 3 include meson condensates plus ordinary matter; 5 
include quarks plus other matter (PCL2 also has hyperons). 
r2 < 3.5 and < 2.5 for all EOSs with hyperons, meson 
condensates, and/or quark cores. The shaded region corre- 
sponds to incompatible values of pi and Fi, as discussed in 
the text. 



straints from the simultaneous measurement of mass and 
moment of inertia and of mass and radius. We exhibit in 
subsection F the combined constraint imposed by causal- 



ity, maximum observed mass, and a future moment-of- 
inertia measurement of a star with known mass. 

In exhibiting the constraints, we show a region of the 
4-dimensional parameter space larger than that allowed 
by the presumed uncertainty in the EOS - large enough, 
in particular, to encompass the 34 candidate EOSs con- 
sidered above. The graphs in Fig. |4] display the ranges 
lO^^ Myne/cm^ < Pi < 103^■^dyne/cm^ 1.4 < Fi < 5.0, 
1.0 < F2 < 5.0, and 1.0 < F3 < 5.0. Also shown is the 
location in parameter space of each candidate EOS, de- 
fined as the set of parameters that provide the best fit 
to that EOS. The shaded region in the top graph corre- 
sponds to incompatible values of pi and Fi mentioned in 
Sect. |IV| These are values for which the pressure pi is so 
large and Fi so small that no curve p{p) can start from 
the low-density EOS above neutron drip and reach pi at 
pi unless the slope ci log p/d log p exceeds Fi. 

To find the constraints on the parameterized EOS im- 
posed by the maximum observed mass and spin, one 
finds the maximum mass and spin of stable neutron stars 
based on the EOS associated with each point of param- 
eter space. A subtlety in determining these maximum 
values arises from a break in the sequence of stable equi- 
libria " an island of unstable configurations - for some 
EOSs. The unstable island is typically associated with 
phase transitions in a way we now describe. 

Spherical Newtonian stars described by EOSs of the 
form p = p(p) are unstable when an average value F of 
the adiabatic index falls below 4/3. The stronger-than- 
Newtonian gravity of relativistic stars means that insta- 
bility sets in for larger values of F, and it is ordinarily this 
increasing strength of gravity that sets an upper limit on 
neutron-star mass. EOSs with phase transitions, how- 
ever, temporarily soften above the critical density and 
then stiffen again at higher densities. As a result, con- 
figurations whose inner core has density just above the 
critical density can be unstable, while configurations with 
greater central density can again be stable. Models with 
this behavior are considered, for example, by Glenden- 
ning and Kettner Bejger et al. jH] and by Zdunik et 
al. [7J (these latter authors, in fact, use piecewise poly- 
tropic EOSs to model phase transitions). 

For our parameterized EOS, instability islands of this 
kind can occur for F2 < 2, when Fi > 2 and F3 > 2. A 
slice of the four-dimensional parameter space with con- 
stant Fi and F3 is displayed in Fig.|5] The shaded region 
corresponds to EOSs with islands of instability. Con- 
tours are also shown for which the maximum mass for 
each EOS has the constant value 1.7A/0 (lower contour) 
and 2.OM0 (upper contour). 

An instability point along a sequence of stellar models 
with constant angular momentum occurs when the mass 
is maximum. On a mass-radius curve, stability is lost in 
the direction for which the curve turns counterclockwise 
about the maximum mass, regained when it turns clock- 
wise. In the bottom graph of Fig. |5] mass-radius curves 
are plotted for six EOSs, labeled A-F, associated with 
six correspondingly labeled EOSs in the top figure. The 
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sequences associated with EOSs B, C and E have two 
maximum masses (marked by black dots in the lower fig- 
ure) separated by a minimum mass. As one moves along 
the sequence from larger to smaller radius - from lower 
to higher density, stability is temporarily lost at the first 
maximum mass, regained at the minimum mass, and per- 
manently lost at the second maximum mass. 

It is clear from each graph in Fig. |5] that either of 
the two local maxima of mass can be the global maxi- 
mum. On the lower boundary (containing EOSs A and 
D), the lower density maximum mass first appears, but 
the upper-density maximum remains the global maxi- 
mum in a neighborhood of the boundary. Above the 
upper boundary (containing EOS F), the higher-density 
maximum has disappeared, and near the upper boundary 
the lower-density maximum is the global maximum. 



A. Causality 

For an EOS to be considered physically reasonable, 
the adiabatic speed of sound cannot exceed the speed of 
light. Perfect fluids have causal time evolutions (satisfy 
hyperbolic equations with characteristics within the light 
cone) only if Vs (the phase velocity of sound) is less than 
the speed of light. An EOS is ruled out by causality if 
Ds > 1 for densities below the central density pmax of 
the maximum-mass neutron star for that EOS. An EOS 
that becomes acausal beyond Pmax at density higher than 
this can always be altered for p > Pmax to a causal EOS. 
Because the original and altered EOS yield identical se- 
quences of neutron stars, causality should not be used 
to rule out parameters that give formally acausal EOSs 
above pmax- 

We exhibit the causality constraint in two ways, first by 
simply requiring that each piecewise polytrope be causal 
at all densities and then by requiring only that it be 
causal below Pmax- The first, unphysically strong, con- 
straint, shown in Fig. [6] is useful for an intuitive under- 
standing of the constraint: The speed of sound is a mea- 
sure of the stiffness of the EOS, and requiring causality 
eliminates the largest values of F^ and pi . 

Fig. [7] shows the result of restricting the constraint to 
densities below Pmax, with the speed of sound given by 
Eq. pO| . A second surface is shown to account for the 
inaccuracy with which a piecewise polytropic approxima- 
tion to an EOS represents the speed of sound. In all but 
one case (BGNIHI) the fits to the candidate EOSs over- 
predict the maximum speed of sound, but none of the 
fits to the candidate EOSs mispredict whether the can- 
didate EOS is causal or acausal by more than 11% (frac- 
tional difference between fit and candidate). We adopt 
as a suitable causality constraint a restriction to a region 
bounded by the surface Ws,max = l-l-mean-l-lcr — 1.12, 
corresponding to the mean plus one standard deviation 
in the error between Vs max for the candidate and best fit 
EOSs. 

In the lower parts of each graph in Fig. [7] where 
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FIG. 5: The region in parameter space where two stable 
neutron-star sequences can occur is shaded in the top figure. 
Contours of constant maximum mass are also shown. The 
higher central density maximum mass contour is solid while 
the lower central density maximum mass contour is dashed. 
Mass-radius curves are plotted for several EOSs in the bottom 
figure. Although difficult to see, EOS C does in fact have a 
second stable sequence. 



Pi < 10^^ dyne/cm^, the bounding surface has the char- 
acter of the first causality constraint, with the restriction 
on each of the three variables pi , r2 and F3 becoming 
more stringent as the other parameters increase, and with 
F3 restricted to be less than about 3. In this low-pressure 
part of each graph, the surface is almost completely inde- 
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pendent of the value of Fi: Because the constraint takes 
the form Tip/{e + p) < (for p <C e) and p < pi is 
so low, the constraint rules out values of Fi only at or 
beyond the maximum Fi we consider. 

In the upper part of each graph, where pi > 
10"^^ dyne/cm^, unexpected features arise from the fact 
that we impose the causality constraint only below the 
maximum density of stable neutron stars - below the 
central density of the maximum-mass star. 

The most striking feature is the way the constraint sur- 
face turns over in the upper part of the top graph, where 
Pi > 10'^^ dyne/cm^, in a way that allows arbitrarily 
large values ofpi. This occurs because, when pi is large, 
the density of the maximum-mass star is small, and a vi- 
olation of causality typically requires high density. That 
is, when the density is low, the ratio p/(e + p) in Eq. ( 10 ) 
is small. As a result, in the top graph, Vs remains too 
small to violate causality before the maximum density is 
reached. In the bottom graph, with Fi — 3.8, Fi is now 
large enough in Eq. (10 1 that the EOS becomes acausal 
just below the transition to F2. This is the same effect 
that places the upper limit on pi seen in the second graph 
of Fig. [6] 

A second feature of the upper parts of each graph is 
the exact independence of the bounding surface on F3. 
The reason is simply that in this part of the parameter 
space the central density of the maximum mass star is 
below p2, implying that no stable neutron stars see F3. 

Finally, we note that in both graphs, for small F2 (the 
right of the graph) , the EOSs yield the sequences men- 
tioned above, in which an island of instability separates 
two stable sequences, each ending at a local maximum of 
the mass. Requiring fs,max to satisfy causality for both 
stable regions rules out EOSs below the lower part of the 
bifurcated surface. 



ri=2.i 




FIG. 6: Causality constraints are shown for two values of Fi. 
For each EOS in the parameter space the maximum speed of 
sound over all densities is used. The shaded surface separates 
the EOS parameter space into a region behind the surface 
allowed by causality (labeled Ws.max < 1) and a region in which 
corresponding EOSs violate causality at any density (labeled 

^s,max ^ 1). 



B. Maximum Mass 



A stringent observational constraint on the EOS pa- 
rameter space is set by the largest observed neutron- 
star mass. Unfortunately, the highest claimed masses 
are also subject to the highest uncertainties and system- 
atic errors. The most reliable measurements come from 
observations of radio pulsars in binaries with neutron 
star companions. The masses with tightest error bars 
(about 0.01 Mq) cluster about 1.4 Mq [Mj. Recent ob- 
servations of millisecond pulsars in globular clusters with 
non-neutron star companions have yielded higher masses: 
Ter 51 and Ter 5 J (35], MSB [36^, PSR J1903-f0327 [37], 
and PSR J0437-4715 [31] all have 95% confidence lim- 
its of about 1.7 Mq, and the corresponding limit for 
NGC 6440B [31] is about 2.3 Mq. However these sys- 
tems are more prone to systematic errors: The pulsar 
mass is obtained by assuming that the periastron ad- 
vance of the orbit is due to general relativity. Perias- 
tron advance can also arise from rotational deformation 
of the companion, which is negligible for a neutron star 



but could be much greater for pulsars which have white 
dwarf or main sequence star companions. Also the mass 
measurement is affected by inclination angle, which is 
known only for the very nearby PSR J0437-4715. And 
with the accumulation of observations of these eccentric 
binary systems (now about a dozen) it becomes more 
likely that the anomalously high figure for NGC 6440B 
is a statistical fluke. Fig. [8] shows the constraint on the 
EOS placed by the existence of 1.7 Mq neutron stars, 
which we regard as secure. Also shown in the figure are 
the surfaces associated with maximum masses of 2.0 Mq 
and 2.3 Mq. 

Since all of the candidate high-mass pulsars are spin- 
ning slowly enough that the rotational contribution to 
their structure is negligible, the constraint associated 
with their observed masses can be obtained by computing 
the maximum mass of nonrotating neutron stars. Cor- 
responding to each point in the parameter space is a se- 
quence of neutron stars based on the associated parame- 
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FIG. 7: Causality constraint as in Fig. [6] However, here, only 
the maximum speed of sound up to the central density of the 
maximum mass star is considered. A second, outlined surface 
shows a weaker constraint to accommodate the expected error 
in the speed of sound associated with a piecewise polytropic 
approximation to an EOS. With a the standard deviation in 
between an EOS and its parameterized representation, 
as measured by the collection of candidate EOSs, the outlined 
surface depicts Vs.max = l+mean+la = 1.12 constraint. 



terized EOS; and a point of parameter space is ruled out 
if the corresponding sequence has maximum mass below 
the largest observed mass. We exhibit here the division 
of parameter space into regions allowed and forbidden by 
given values of the largest observed mass. 

We plot contours of constant maximum mass in Fig. [8] 
Because EOSs below a maximum mass contour produce 
stars with lower maximum masses, the parameter space 
below these surfaces is ruled out. The error in the maxi- 
mum mass between the candidate and best fit piecewise 
polytropic EOSs is |mean|-|-l(T = 1.7% (the magnitude of 
the mean error plus one standard deviation in the error 
over the 34 candidate EOSs), so the parameters that best 
fit the true EOS are unlikely to be below this surface. 

The surfaces of Fig. [8] have minimal dependence on Fi, 
indicating that the maximum mass is determined primar- 
ily by features of the EOS above pi. In Fig. [8] we have set 



Fi to the least constraining value in the range we consider 
- to the value that gives the largest maximum mass at 
each point in {pi,F2,F3} space. Varying Fi causes the 
contours to shift up, constraining the parameter space 
further, by a maximum of 10*^'^ dyne/cm^. The depen- 
dence of the contour on Fi is most significant for large 
values of pi where the average density of a star is lower. 
The dependence on Fi decreases significantly as pi de- 
creases. 




FIG. 8; The above surfaces represent the set of parameters 
that result in a constant maximum mass. An observation 
of a massive neutron star constrains the equation of state to 
lie above the corresponding surface. Fi is set to the least 
constraining value at each point. The lower shaded surface 
represents Mmax = 1.7 Mq; the middle and upper (outlined) 
surfaces represent Mmax = 2.0 Mq and Mmax ~ 2.3 Mq 
respectively. 



As discussed above, some of the EOSs produce se- 
quences of spherical neutron stars with an island of insta- 
bility separating two stable sequences, each with a local 
maximum of the mass. As shown in Fig. [5] this causes a 
contour in parameter space of constant maximum mass 
to split into two surfaces, one surface of parameters which 
has this maximum mass at the lower pc local maximum 
and another surface of parameters which has this maxi- 
mum mass at higher pc branches. Since such EOSs allow 
stable models up to the largest of their local maxima, 
we use the least constraining surface (representing the 
global maximum mass) when ruling out points in param- 
eter space. 



C. Gravitational redshift 

We turn next to the constraint set by an observed 
redshift of spectral lines from the surface of a neutron 
star. We consider here only stars for which the broad- 
ening due to rotation is negligible and restrict our dis- 
cussion to spherical models. The redshift is then z = 
(1 - 2M/i?)"i/2 _ 1^ and measuring it is equivalent to 
measuring the ratio M/R. With no independent mea- 
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surement of mass or radius, the associated constraint 
again restricts the parameter space to one side of a sur- 
face, to the EOSs that allow a redshift as large as the 
largest observed shift. For spherical models, the con- 
figuration with maximum redshift for a given EOS is or- 
dinarily the maximum-mass star. By increasing pi, T2 or 
Fs, one stiffens the core, increasing the maximum mass, 
but also increasing the radius at fixed mass. The out- 
come of the competition usually, but not always, yields 
increased redshift for larger values of these three param- 
eters; that is, the increased maximum mass dominates 
the effect of increased radius at fixed mass for all but the 
largest values of pi . 

Cottam, Paerels, and Mendez [H] claim to have ob- 
served spectral lines from EXO 0748-676 with a gravi- 
tational redshift of z = 0.35. With three spectral lines 
agreeing on the redshift, the identification of the spec- 
tral features with iron lines is better founded than other 
claims involving only a single line. The identification re- 
mains in doubt, however, because the claimed lines have 
not been seen in subsequent bursts [15]. There is also a 
claim of a simultaneous mass-radius measurement of this 
system using Eddington-limited photospheric expansion 
x-ray bursts |43j which would rule out many EOSs. This 
claim is controversial, because the 95% confidence inter- 
val is too wide to rule out much of the parameter space, 
and we believe the potential for systematic error is under- 
stated. However, the gravitational redshift is consistent 
with the earlier claim of 0.35. Thus we treat z — 0.35 as a 
tentative constraint. We also exhibit the constraint that 
would be associated with a measurement of z = 0.45. 



Our parameterization can reproduce the maximum 
redshift of tabulated EOSs to 3.2% (mean-|-lcr). Fig- 
ure [9] displays surfaces of constant redshift z = 0.35 and 
2: = 0.45 for the least constraining value of Fi = 5 in the 
range we consider. Surfaces with different values of Fi are 
virtually identical for pi < lO'^^-^ dyne/cm^, but diverge 
for higher pressures when Fi is small (< 2.5). In the dis- 
played parameter space, points in front of the z = 0.35 
surface, corresponding to stiffer EOSs in the inner core, 
are allowed by the potential z = 0.35 measurement. From 
the location of the z = 0.35 and z = 0.45 surfaces, it is 
clear that, without an upper limit on Fi < 2.5, an ob- 
served redshift significantly higher than 0.35 is needed to 
constrain the parameter space. In particular, most of the 
parameter space ruled out by z = 0.35 is already ruled 
out by the My^ax = 1.7 Mq constraint displayed in Fig.|8] 



^ One could also imagine a measured redshift small enough to 
rule out a class of EOSs. The minimum redshift for each EOS, 
however, occurs for a star whose central density is below nuclear 
density. Its value, 2 5 X 10~^, thus depends only on the EOS 
below nuclear density. (See, for example Haensel et al.|40|.) 




FIG. 9: Surfaces in the EOS parameter space for which the 
maximum redshift of stable spherical neutron stars has the 
values 0.35 (shaded surface) and 0.45 (outlined surface). A 
measured redshift from the surface of a neutron star would 
exclude the region of parameter space behind the correspond- 
ing surface. Fi is fixed at 5.0, the least constraining in the 
range we considered. 



D. Maximum Spin 

Observations of rapidly rotating neutron stars can also 
constrain the EOS. The highest uncontroversial spin fre- 
quency is observed in pulsar Ter 5AD at 716 Hz [33]. 
There is a claim of 1122 Hz inferred from oscillations in 
x-ray bursts from XTE J 1239-285 [45], but this is con- 
troversial because the statistical significance is relatively 
low, the signal could be contaminated by the details of 
the burst mechanism such as fallback of burning material, 
and the observation has not been repeated. 

The maximum angular velocity of a uniformly rotat- 
ing star occurs at the Kepler or mass-shedding hmit, Qk, 
with the star rotating at the speed of a satellite in circular 
orbit at the equator. For a given EOS, the configuration 
with maximum spin is the stable configuration with high- 
est central density along the sequence of stars rotating at 
their Kepler limit. An EOS thus maximizes rotation if 
it maximizes the gravitational force at the equator of a 
rotating star - if it allows stars of large mass and small 
radius. To allow high mass stars, the EOS must be stiff 
at high density, and for the radius of the high-mass con- 
figuration to be small, the EOS must be softer at low 
density, allowing greater compression in the outer part 
of the star [321 137]. In our parameter space, a high an- 
gular velocity then restricts one to a region with large 
values of F2 and F3, and small values of pi and Fi. 

As with the maximum mass, the maximum frequency 
is most sensitive to the parameter pi, but the frequency 
constraint complements the maximum mass constraint 
by placing an upper limit on pi over the parameter space, 
rather than a lower limit. 

To calculate the maximum rotation frequencies for our 
parameterized EOS, we used the open-source code rns 
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for axisymmetric rapid rotation in the updated form 
rns2.0 |48]. For a given EOS, the model with maxi- 
mum spin is ordinarily close to the model with maximum 
mass, but that need not be true for EOSs that yield two 
local mass maxima. The resulting calculation of maxi- 
mum rotation requires some care, and the method we use 
is described in Appendix B. The error incurred in using 
the parameterized EOS instead of a particular model is 
2.7% (mean+lcr). 

Spin frequencies of 716 Hz and even the possible 
1122 Hz turn out to be very weak constraints because 
both are well below the Kepler frequencies of most EOSs. 
Thus we plot surfaces of parameters giving maximum 
rotation frequencies of 716 Hz in Fig. [TO] and 1300 Hz 



and 1500 Hz in Fig. 11 The region of parameter space 
above the maximum observed spin surface is excluded. 
In the top figure, maximum mass stars have central den- 
sities below p2 so there is no dependence on F3. In the 
bottom figure the least constraining value of Fi = 5 is 
fixed. The surface corresponding to a rotation of 716 
Hz only constrains the parameter space that we consider 
(pi < 10^^-^ dyne/cm^) if Fi < 2.5. The minimum 
observed rotation rate necessary to place a firm upper 
limit on pi is roughly 1200 Hz for Fi = 5. The surface 
/max = 1500 Hz for Fi = 5 is also displayed in Fig.[TT]to 
demonstrate that much higher rotation frequencies must 
be observed in order to place strong limits on the param- 
eter space. 




FIG. 10: The above surface represents the set of parameters 
that result in a maximum spin frequency of 716 Hz for the top 
surface. For high values of pi there is no dependence on T3. 
The wedge at the back riglit is the shaded region of Fig. |4] 
corresponding to incompatible values of pi and Fi. 



Because it is computationally expensive to use rns 
to evaluate the maximum rotation frequency for a wide 
range of values in a 4-parameter space, one can also use 
an empirical formula. Haensel and Zdunik |49| found 
that the maximum stable rotation for a given EOS can 
be found from the maximum-mass spherically symmetric 




FIG. 11: The above surfaces represent the set of parameters 
that result in a maximum spin frequency of 1200 Hz for the 
top surface and 1500 Hz for the bottom surface. That is, 
observations of such high spin frequencies would constrain 
the EOS to lie below the corresponding surface. For these 
surfaces Fi = 5, the least constraining value. 



model for that EOS with mass il/g and radius i?s 

/ ^''max \ / -''^^s \ / -'^s 



Vio^s- 



Mq 



10 km 



(13) 



In other words the maximum rotation is proportional to 
the square root of the average density of the star. 

The original calculation of Haensel and Zdunik gave 
K = 0.77. An overview of subsequent calculations is 
given by Haensel et al. in |50j . reporting values of 
K ~ 0.76 — 0.79 for a range of EOS sets and calcula- 
tion methods including those of [5TJ [5H [53] . If we cal- 
culate maximum rotations with rns as described above, 
using the 34 tabulated EOSs, we find k = 0.786 ± 0.030. 
The corresponding best fit parameterized EOSs give k — 
0.779 ± 0.027. 



E. Moment of inertia or radius of a neutron star of 
known mass 

The moment of inertia of the more massive component, 
pulsar A, in the double pulsar PSR J0737-3039 may be 
determined to an accuracy of 10% within the next few 
years [T^ by measuring the advance of the system's pe- 
riastron, and implications for candidate EOSs have been 
examined in [H1|S31[5S] . As noted earlier, by finding both 
mass and moment of inertia of the same star one imposes 
a significantly stronger constraint on the EOS parameter 
space than the constraints associated with measurements 
of mass or spin alone: The latter restrict the EOS to the 
region of parameter space lying on one side of a surface, 
the region associated with the inequality Mi„ax(pi,ri) > 

-^observed Or with ^ma.x{Ply^i) > ^^observed- The simul- 
taneous measurement, on the other hand, restricts the 
EOS to a single surface. That is, in an n-dimensional pa- 
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rameter space, the full n-dimensional set of EOSs which 
allow a 1.338 Mq model, and those EOSs for which 
that model has moment of inertia /observed form the 
{n — l)-dimensional surface in parameter space given by 
I{pi,Ti,M — 1.338Af0) = /observed- (Wc usc here the 
fact that the 44 Hz spin frequency of pulsar A is slow 
enough that the moment of inertia is nearly that of the 
spherical star.) Moreover, for almost all EOSs in the pa- 
rameter space, the central density of a 1.338 Mq star 
is below the transition density p^- Thus the surfaces of 
constant moment of inertia have negligible dependence 
on Fa, the adiabatic index above p2, and the EOS is re- 
stricted to the two-dimensional surface in the pi-ri-r2 
space given by /(pi, Ti, M = I.338M0) = /observed- 

This difference in dimensionality means that, in prin- 
ciple, the simultaneous equalities that give the constraint 
from observing two features of the same star are dramati- 
cally stronger than the inequalities associated with mea- 
surements of mass or spin alone. In practice, however, 
the two-dimensional constraint surface is thickened by 
the error of the measurement. The additional thickness 
associated with the error with which the parametrized 
EOS can reproduce the moment of inertia of the true 
EOS is smaller, because the parameterized EOS repro- 
duces the moment of inertia of the 34 candidate EOSs to 
within 2.8% (|mean| -I- la). 

In Fig. [12] we plot surfaces of constant moment of in- 
ertia that span the range associated with the collection 
of candidate EOSs. The lower shaded surface represents 
/ = 1.0 X 10**^ gcm^. This surface has very little depen- 
dence on Fi because it represents a more compact star, 
and thus for a fixed mass, most of the mass is in a denser 
state p > pi. The structures of these stars do depend on 
F3, and the corresponding dependence of / on F3 is shown 
by the separation between the surfaces in Fig. [12] The 
middle outlined surface represents / = 1.5 x 10*'' gcm^, 
and is almost a surface of constant pi . The top outlined 
surface represents / = 2.0 x 10^^ gcm^. This surface has 
little dependence on F2, because a star with an EOS on 
this surface would be less compact and thus most of its 
mass would be in a lower density state p < pi. 

If the mass of a neutron star is already known, 
a measurement of the radius constrains the EOS to 
a surface of constant mass and radius, R{pi,Ti) = 
/^observed, M{pi,Ti) = Mobserved in the 4-dimensional pa- 
rameter space. The thickness of the surface is dominated 
by the uncertainty in the radius and mass measurements, 
since our parameterization produces the same radius as 
the candidate EOSs to within 1.7% (|mean| -t- la). We 
plot in Fig 
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star that span the range of radii associated with the col- 
lection of candidate EOSs. As with the moment of in- 
ertia, the radius depends negligibly on F3 as long as the 
radius is greater than 11 km. For smaller radii, the vari- 
ation with F3 is shown by the separation between the 
surfaces in Fig. [13] 

Very recently analyses of time-resolved spectroscopic 
data during thermonuclear bursts from two neutron stars 




FIG. 12: The above surfaces represent the set of parame- 
ters that result in a star with a mass of 1.338 Mq and a 
fixed moment of inertia, i.e. possible near-future measure- 
ments of PSR J0737-3039A. / = 1.0 x 10*^ gcm^ for the 
shaded surfaces, whose separation corresponds to varying 
Fa. / = 1.5 X 10^^ gcm^ for the middle outlined surface. 
I = 2.0 X 10*^ gcm^ for the top outlined surface. The wedge 
at the back right is the shaded region of Fig. [4] corresponding 
to incompatible values of pi and Fi. 



in low-mass x-ray binaries were combined with distance 
estimates to yield M = 1.4 Mq and i? = 11 km or M = 

1.7 M© and /? = 9 km for EXO 1745-248 [56 and M = 

1.8 Mq and /? = 10 km for 4U 1608-52 [57J,'both with 
error bars of about 1 km in R. These results are more 
model dependent than the eventual measurement of the 
moment of inertia of PSR J0737-6069A, but the accuracy 
of the measurement of / remains to be seen. 




FIG. 13: The above surfaces represent the set of parame- 
ters that result in a star with a mass of 1.4 Mq and a fixed 
radius. R — 9 km for the shaded surfaces, whose separation 
corresponds to varying F3. i? = 12 km for the middle outlined 
surface, i? = 16 km for the top outlined surface. The wedge 
at the back right is the shaded region of Fig. [4] corresponding 
to incompatible values of pi and Fi. 



14 



F. Combining constraints 

The simultaneous constraints imposed by causality, a 
maximum observed mass of I.TMq, and a future mea- 
surement of the moment of inertia of PSR J0737-3039A, 
restrict the parameter space to the intersection of the al- 
lowed regions of Figs. [7] |8] and [12] We show in Fig. [14] 
the projection of this jointly constrained region on the 
P1 — T2— Fa subspace. This allows one to see the cutoffs 
imposed by causality that eliminate large values of F2 
and F3 and (in the top figure) the cutoffs imposed by the 
existence of a I.TMq model that eliminates small values 
of F2 and F3. 

We noted above that measuring the moment of inertia 
of a I.338M0 star restricts the EOS at densities below p2 
to a two-dimensional surface in the pi — Fi — F2 space. In 
the full 4-dimensional parameter space, the correspond- 
ing surfaces of constant M and I of Fig. [l4] are then 
three dimensional and independent of F3. Their projec- 
tions onto the pi — F2 — F3 subspace are again three- 
dimensional and independent of F3, their thickness due 
to the unseen dependence of the mass and moment of 
inertia on Fi. For small moments of inertia there is neg- 
ligible dependence on Fi so the allowed volume in Fig. [14] 
is thin. The thickness of the allowed volume increases as 
the moment of inertia increases because the dependence 
on Fi also increases. 

In Fig. [15] we explore a relation between the moment 
of inertia /(1. 338) of PSR J0737-3039A and the max- 
imum neutron star mass, in spite of the fact that the 
maximum mass is significantly greater than 1.338 Mq. 
For three values of the moment of inertia that span 
the full range associated with our collection of candi- 
date EOSs, we show joint constraints on r2 and F3 in- 
cluding causality and maximum neutron star mass. For 
/(1. 338) = 1.0 X 10**^ gcm^, F2 is nearly unconstrained, 
while F3 is required to lie in a small range between the 
causality constraint and the reliable observations of stars 
with mass 1.7 Mq. For larger values of /(1. 338), F2 is 
more constrained and F3 is less constrained. However, 
the highest values of /(1. 338) are associated with the 
highest maximum neutron star masses. Thus, if a neu- 
tron star mass of about 2.3 AIq is confirmed, it implies 
that /(1.338) is about 2 x 10'*^ gcm^. Conversely iff 
7(1.338) is measured first and is about 1 x 10'''^ gcm^, it 
implies that the maximum neutron star mass is less than 
about 1.9 Mq. 

The allowed range for pi as a function of the moment 
of inertia of J0737-3039A is shown in Fig. [16] The entire 
shaded range is allowed for a 1.7 Mq maximum mass. 
The medium and darker shades are allowed for a 2.0 Mq 
maximum mass. Only the range with the darker shade is 
allowed if a 2.3 Mq star is confirmed. It should be noted 
that for small moments of inertia, this plot overstates the 
uncertainty in the allowed parameter range. As shown in 
Fig. [14] the allowed volume in F2 — F3 — pi space for 
a small moment of inertia observation is essentially two 
dimensional. If the moment of inertia is measured to be 




FIG. 14; The figure portrays the joint constraint imposed by 
causality (i's,max < l+mean-f-lo-), the existence of a 1.7 Mq 
neutron star, and by a future measurement of the moment 
of inertia I of J0737-3039A. Each thick shaded surface is the 
volume in r2 — Fs — pi space allowed by the joint constraint 
for the labeled value of I. 



this small, then the EOS would be better parameterized 
with the linear combination Q;log(pi) -I- PT2 instead of 
two separate parameters log(pi) and F2. 



VII. DISCUSSION 

We have shown how one can use a parameterized piece- 
wise polytropic EOS to systematize the study of observa- 
tional constraints on the EOS of cold, high-density mat- 
ter. We think that our choice of a 4-parameter EOS 
strikes an appropriate balance between the accuracy of 
approximation that a larger number of parameters would 
provide and the number of observational parameters that 
have been measured or are likely to be measured in the 
next several years. The simple choice of a piecewise poly- 
trope, with discontinuities in the polytropic index, leads 
to suitable accuracy in approximating global features of 
a star. But the discontinuity reduces the expected accu- 
racy with which the parameterized EOS can approximate 
the local speed of sound. One can largely overcome the 
problem by using a minor modification of the parameter- 
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FIG. 15; The allowed values of r2 and Fa depend strongly 
on the moment of inertia of PSR J0737-3039A. In top, 
middle and bottom figures, respectively, / has the values 
1.0x10*^ gcm^ I = 1.5x10'*'^ gcm^ and 7 = 2.0x10^^ gcm^ 
In each figure the upper curves are the Us,max = 1 (dotted) and 
fs,max ~ l+mean+l(j — 1.12 (solid) causality constraints. 
Shading indicates a range of possible maximum mass con- 
straints, with increasing maximum mass leading to a smaller 
allowed area. All shaded areas are allowed for a 1.7 Mq max- 
imum neutron star mass. The medium and dark shades are 
allowed if a 2.0 Mq star is confirmed. Only the darkest shade 
is allowed if a 2.3 Mq star is confirmed. 



ized EOS in which a fixed smoothing function near each 
dividing density is used to join the two polytropes. 




1.0 1.2 1.4 1.6 1.8 2.0 
Moment of Inertia (10^*^ g cm^) 



FIG. 16: The allowed range of pi as a function of the mo- 
ment of inertia of J0737-3039A when combined with causality 
(fs,max ~ l-fmean-l-la) and observed mass constraints. All 
shaded areas are allowed by a 1.7 Mq maximum mass. The 
medium and dark shades are allowed if a 2.0 Mq star is con- 
firmed. Only the darkest shade is allowed if a 2.3 AIq star is 
confirmed. 



We see that high-mass neutron stars are likely to pro- 
vide the strongest constraints from a single measure- 
ment. The work dramatizes the significantly more strin- 
gent constraints associated with measurements like this, 
if two (or more) physical features of the same star can be 
measured, and an n-dimensional parameter space is re- 
duced by one (or more) dimension(s), to within the error 
of measurement. In particular, a moment of inertia mea- 
suremement for PSR J0737-3039 (whose mass is already 
precisely known) could strongly constrain the maximum 
neutron star mass. 

The effect of EOS-dependent tidal deformation can 
modify the gravitational waves produced by inspiraling 
neutron stars. This modification is largely dependent 
on the radius of the neutron star. Flanagan and Hin- 
derer [SS] investigate constraints on an EOS-dependent 
tidal parameter, the Love number, from observations 
of early inspiral. A companion to this paper |llj uses 
the parametrized EOS in numerical simulations to ex- 
amine the future constraint associated with expected 
gravitational- wave observations of late inspiral in binary 
neutron stars. 

Finally, we note that the constraints from observations 
of different neutron star populations constrain different 
density regions of the EOS. For moderate mass stars such 
as those found in binary pulsar systems, the EOS above 
P2 = 10^^ '' g/cm^ is unimportant. For near-maximum 
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mass stars, the EOS below pi = 10^^'^ g/cm^ has little 
effect on neutron star properties. This general behavior 
is independent of the details of our parameterization. 
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APPENDIX A: EVALUATING MASS, RADIUS, 
AND MOMENT OF INERTIA 



The moment of inertia of a rotating star is the ra- 
tio / = J/f2, with J the asymptotically defined angular 
momentum. In finding the moment of inertia of spheri- 
cal models, we use Hartle 's slow-rotation equations [55] . 
adapted to piecewise polytropes in a way we describe be- 
low. The metric of a slowly rotating star has to order Q. 
the form 

+r^d6'^ + r^sm^ed(j)'^, (Al) 

where and A are the metric functions of the spherical 
star, given by 
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The frame-dragging uj{r) is obtained from the t(f> compo- 
nent of the Einstein equation in the form 



)^d_ 
r^ dr 



i.du}\ 
''Tr) 



Ado_ 
r dr 



= 0, 



(A6) 



where w = f2 — w is the angular velocity of the star mea- 
sured by a zero-angular-momentum observer and 
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(A7) 



The angular momentum is obtained from lo, which has 
outside the star the form lo = 2J /r^ . 

In adapting these equations, we roughly follow Lind- 
blom [60j, replacing r as a radial variable by a general- 
ization rj := h — \ oi the Newtonian enthalpy. ^ Because 
77 is monotonic in r, one can integrate outward from its 
central value to the surface, where 77 = 0. 

For the piecewise polytropes of Sec. |III| the equation 
of state given in terms of 77 is 
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(A9) 
(AlO) 



where = l/(ri — 1) is the polytropic index. 

This replacement exploits the first integral /le* = 
^/l — 2M/R of the equation of hydrostatic equilibrium 
to eliminate the differential equation (A3 ) for <i>; and the 



enthalpy, unlike e and p, is smooth at the surface for a 
polytropic EOS. Eqs. ( A4||A6 1 are then equivalent to the 
first-order set 
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where a dZu/dr. 

The integration to find the mass, radius, and moment 
of inertia for a star with given central value rj — rjc 
proceeds as follows: Use the initial conditions r{rjc) = 
"^{ic) = OL{rjc) = and arbitrarily choose a central value 
LjQ of Lj. Integrate to the surface where r] — Q, to obtain 
the radius R = r{r] = 0) and mass M = m^rj = 0). The 
angular momentum J is found from the radial derivative 
of the equation 



2J 



(A15) 



* Lindblom, however, uses log h instead of — 1 as his radial vari- 
able. Because of the form of the piecewise polytrope, h — 1 is 
a more convenient choice here. This variable is also used by 
Haensel and Potekhin in [9]. 
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The mass-shed hmit gives a maximally rotating equi- 
librium model for each central energy density Cc, but, as 
in the spherical case, these equilibrium models are not 
guaranteed to be stable to perturbations. 

Overall stability in uniformly rotating models is gov- 
erned by the stability of the model to pseudoradial per- 
turbations. As in the spherical = case, there can 
exist alternating regions of stable and unstable rotating 
models along a sequence of fixed fl. A criterion for the 
onset of instability is developed by Friedman, Ipser and 
Sorkin in |61| : The critical points that potentially indi- 
cate a change in stability are extrema of mass-energy M 
under variation in both baryon mass Aff, and angular mo- 
mentum J; and can be determined by extremizing Mh on 
sequences of constant J or extremising J on sequences 
of constant Mi,. Universally valid searches for limiting 
stability, as in for example have therefore required 
explicitly covering the set of models with sequences of 
constant rest mass and extremizing J on each one, 
or vice versa — a computationally expensive procedure. 

For most npe/i-only EOSs, the maximally rotating sta- 
ble model is close to the point on the mass-shed limit with 
maximal mass-energy, and this model has been used for 
an estimate of maximal rotation in surveys of large num- 
bers of EOSs. However, this is not always the case. An 
example is in EOS L of [62], or the parameterized EOS 
of Fig. [17| 

Consider the two-parameter family of rotating neutron 
stars as a surface S in M^-J-ec space. The central energy 
density tc and axis ratio r are suitable parameters for this 
surface. At points where Mf, is maximum along constant 
J sequences, the vector tangent to the sequence points in 
the Cc direction. Models with limiting stability are found 
where the tangent plane to the surface of equilibrium 
models contains a vector in the Cc direction, e^. 

Given the parameterization of the surface in terms of 
Cc and r, the normal vector to the surface E — {A/fc, J, ed 
is along 



FIG. 17: Family of rotating models for a given EOS, from 
the spherical limit to the mass-shed limit. The surface is pro- 
jected into Mb-ec and covered with lines of constant J. The 
maximum Mi, on each line gives the critically stable model; 
rotation at the mass-shed limit increases as tc increases past 
the maximum mass model. 



with component along q,. 

dMb dJ dJ dMb 

Tl ^= — " 

'^'^ dcr dx de^ dx 



(B2) 



which is zero at the critical line between stable and unsta- 
ble equilibriums on the surface S. A covariant statement 
of this condition for marginal stability is dAff, /\ dJ — 0. 

The maximally rotating model for a given EOS may be 
determined, without finding sequences of constant J and 
Mb, by considering a sequence of central energy densities 
£(,. First, increase the axis ratio r until the Kepler limit 
is found, as in the example program main . c of rnsv2 . 0. 
Second, vary ec and r around this point to estimate the 
partial derivatives of Eq. B2 The sign of n^^ will change 



as 9S 

d€r dt 



(Bl) 



as the Kepler limit sequence crosses the stability limit. 



APPENDIX C: ANALYTIC FITS TO 
TABULATED EOSS 

As another measure of the ability of the parameterized 
EOS to fit candidate EOSs from the literature, we exam- 
ine how well the parameterized EOS reproduces neutron 
star properties predicted by the candidate EOSs. We 
use an analytic form of the (SLy) low-density EOS that 
closely matches its tabulated values. With rms resid- 
ual less than 0.03, p{p) for SLy is approximated between 
p — lO'^ g/cnx^ and p — 10^^ g/cm'^ by four polytropic 
pieces. The four regions correspond roughly to a nonrel- 
ativistic electron g relativistic electron gas, neutron 
drip, and the density range from neutron drip to nuclear 
density. Using the notation of Sect. |III| the analytic form 
of the SLy EOS is set by the values of Ki,Ti and pi listed 



19 



in Table [TT] The parameters for the three piece polytropic 
high-density EOS, the corresponding residuals, as well as 
the observable properties of these EOSs and the error in 
using the best fit parameterized EOS instead of the tab- 



ulated EOS are listed in Table III The parameterized 
EOS systematically overestimates the maximum speed 
of sound. 



TABLE II: An analytic representation of p(p) for the SLy 
EOS below nuclear density uses polytropes specified by the 
constants listed here. Ti is dimensionless, pi is in g/cm"^, and 
Ki is in cgs units for which the corresponding value of p is in 
units of dyne/cm'^. The last dividing density is the density 
where the low density EOS matches the high density EOS 
and depends on the parameters pi and Fi of the high density 
EOS. 



Kj Tj Pi 

6.80110e-09 1.58425 2.44034e-h07 
1.06186e-06 1.28733 3.78358e-hll 
5.32697e4-01 0.62223 2.62780e+12 
3.99874e-08 1.35692 
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TABLE III: Comparison of candidate EOSs and their best fits. The parameters that provide the best fit to the candidate 

EOSs as well as the residual are given, pi is in units of dyne/cm^. Values for observables calculated using the tabulated 
EOSs are also given. Ws,max is the maximum adiabatic speed of sound below the central density of the maximum mass neutron 
star. A'/,nax is the maximum nonrotating mass configuration in units of Mq. Zmax is the corresponding maximum gravitational 
redshift. /max is the maximum rotation frequency in Hz, as calculated using the rotating neutron-star code rns. /1.338 is the 
moment of inertia for a 1.338 Mq star in units of 10*^ gcm^. Ria is the radius of a 1.4 Mq star in units of km. The difference 
in calculated values for each observable when using the tabulated EOS (Otnb) versus the best fit parameterized EOS {Om) is 
calculated with (Ofit/Otab — 1)100. The last entry gives the mean and standard deviation of the errors for each observation. 
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Mean error 

Standard deviation of error 


5.68 
5.52 


-0.71 
0.96 


0.74 
2.42 


0.43 
2.25 


-1.26 
1.57 


-0.37 
1.29 



"The tables for GSl, GS2, and H6 do not go up to the central 
density of the maximum mass star. For most observables, the EOS 
can be safely extrapolated to higher density with minimal error. 
However, the maximum speed of sound is highly sensitive to how 
this extrapolation is done. Thus, we only use the maximum speed 
of sound up to the last tabulated point when comparing the values 
for the table and fit. 

'GSl has a maximum mass less than 1.4 Mq. 



